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We theoretically study classical thermal activation (TA) and macroscopic quantum tun- 
neling (MQT) for a YBa2Cu307_^ (YBCO) Josephson junction coupled with an LC 
circuit. The TA and MQT escape rate are calculated by taking into account the two- 
dimensional nature of the classical and quantum phase dynamics. We find that the MQT 
escape rate is largely suppressed by the coupling to the LC circuit. On the other hand, 
this coupling leads to the slight reduction of the TA escape rate. These results are rel- 
evant for the interpretation of a recent experiment on the MQT and TA phenomena in 
YBCO bi-epitaxial Josephson junctions. 

Keywords: Josephson effects; macroscopic quantum tunneling; high-T c superconductor. 

1. Introduction 



Macroscopic quantum tunneling (MQT) has become a focus of interest in physics 
and chemistry because it can provide a signature of quantum behavior in a macro- 
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scopic system.^ Among several works on MQT, Josephson junctions have been in- 
tensively studied.^ Heretofore, however, experimental and theoretical investigations 
of MQT have been focused on low-T c superconductor Josesphson junctions. 

Renewed interest in MQT occurred after the recent experimental observa- 
tions of MQT in high-T c superconductor Jose phso n junctions, e.g., YBa2Cu307_a 
(YBCO) grain-boundar y bi-epitaxial junctions 3 -^ an( j Bi2Sr 2 CaCu208+,5 intrinsic 
junctions! ^ ' 7 ^ 9 ' ^"^ Such fascinating results open up the possibility for realizing 
high-T c quantum bits (qubits). Over the past year much progress has been achieved 
for maki ng the ory of the MQ ^ 12 ' 13 ' 14 ! 15 ' 16 ' 17 ' 18 ' 19 ] and macroscopic quantum 
coherenc e! 2 ^ 2 ^ in high-T c junctions. 

In YBCO grain boundary junctions which was used in MQT experiments} 3 -^ it 
was found that the stray capacitance Cs of the electrodes, due to the large dielectric 
constant e of the STO substrate at low temperature (e > 10000), and the induc- 
tance L$, due to the large London penetration depth in c-axis and/or Josephson 
coupling between Cu02 planes in one of the electrodes have large influence on the 
macroscopic dynamics and can be taken into account by an extended circuit model 
(see Fig. 1(a)). In Fig. 1, <\> is the phase difference across the Josephson junction 
and 4> s = (2tt/$o)IsLs + 4> is the phase difference across the capacitor Cs, where 
Is is the current through the inductor. As will be mentioned later, the addition 
of the LC circuit results in a two-dimensional potential U(<p, 4>s) which makes the 

dynamics much more complex than for an ordinary single junction (see Fig.l (b)). 

l4l 

It was found that, in the microwave- assisted MQT experiment, - the bias current 
Text depende nce of t he Josephson plasma frequency lo p is quantitatively explained 
by this modelP^2-3J However, the validity of the extended circuit model for thermal 
activation (TA) and MQT escape processes have not yet been explored. In this pa- 
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Fig. 1. (a) The extended circuit model for a bi-epitaxial YBCO junction, including the bias 
current Text, the stray capacitance Cs and kinetic inductance L$- <j> an d <t>S are the phase difference 
across the Josephson junction and the stray capacitance Cs, respectively, (b) The two-dimensional 
potential profile U(<j>, 4>s) °f ( a ) f° r 7 = lext/Ic = 0-92 and rj = Ls/Ljq = 7. 
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per, we investigate the validity of the extended circuit model for both the TA and 
MQT escape processes and discuss the role of the two-dimensiona l natu re on the 
macroscopic quantum dynamics in the YBCO bi-epitaxial junctionP^ESl 

2. Model and Lagrangian 

In this section we derive the Lagrangian for the Josephson junction coupled to the 
LC circuit and discuss the two-dimensional potential structure of this model. The 
Hamiltonian of the circuit [Fig. 1(a)] can be written as 



where Qj = Cj($ /2Tr)(d(j)/dt) and Qs = Cs($o/2-k) {d<ps/dt) are the charge 
on the junction and the stray capacitor Cs, respectively, Ej — (H/2e)Ic is the 
Josephson coupling energy, and <&o = h/2e is the flux quantum. 
The Lagrangian is given by 

Cj /<& o 90\* C S (* Q d<t> s V 



u{4>, <}>s) = Ej 



-cos0+ 70s 

2i] 



(3) 



where 7 = I ext /Ic and rj = 2nI c L s /§ = L S /L JQ with L J0 = $ /2nI c being 
the zero-bias Josephson inductance. The Lagrangian C describes the quantum dy- 
namics of a fictive particle moving in a two-dimensional tilted washboard potential 
£/(0, 0s). Therefore the escape paths traverse a two-dimensional land scape, so the 
macroscopic dynamics in this model becomes more complicated than that in the 
simple one-dimensional model. 

The two-dimensional potential profile t/(0, 0s) is shown in Fig. 1(b). The mean 
slope along the 4>s direction is proportional to the bias current Text, and the mean 
curvature in the diagonal direction (0 = 0s) is due to the inductive coupling between 
the Josephson junction and capacitance Cs characterized by 77. The local minimum 
point (0 m , 0c?) and the saddle point (0*, <fi s ) is given from Eq. ([3]) as 

(4> m ,4>s) = (sin" 1 7, 177 + sin" 1 7), (4) 

0s) = ~ sin ~ 1 7)^ - sin -1 7 + 777) , (5) 
Then the potential barrier height Vq B is given by 

V m = U (0*, S ) - U (0 m , 0S 1 ) = Ej (278m- 1 7 - ^7 + 2 VI - 7*) , (6) 

and is a decreasing function of 7. Importantly, the potential barrier height Vq B does 
not depend on the LC circuit parameters Ls and Cs-^ Moreover, the expression 
of Vq B is identical with the barrier height for usual single Josephson junctions, i.e., 
one-dimensional model.^ 
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Next we will derive an approximate expression of U{4>,(j)s) for 7 < 1. By intro- 



ducing the new coordinate (x, y) = (q 
rewrite the Lagrangian as 

_ M fdx 

T [at 



U(x,y) = Um{x)+Ej 



and assuming 7 < 1, we can 



m 
2 
(x 



dy 
Of 

y? 



2i] 



(7) 
(8) 



where M = Cj(<J> / 27r ) 2 and m = C s {^/2tt) 2 . Here U m (x) « (l/2)Mw ? 2 (i 2 - 
x 3 /xi) is the potential of the Josephson junction without the LC circuit, where 
xi = 3yl — 7 2 and uj p — w p o(l — 7 2 ) 1 ^ 4 is the Josephson plasma frequency with 
u>pQ — y/2irlc /QqC'j being the zero-bias plasma frequency. 



3. Effective Action 

In this section we will derive a n effe ctive action from the Lagrangian ([7]). By using 
the functional integral methodjU^S] ^he partition function Z of the system can be 
written as 



T>x(t) I T>y(r) exp 



drC\x,y\ 



where 



C[x,y] 



-x 2 + —y 2 



U{x,y), 



(9) 



(10) 



2 2 

is the Euclidean Lagrangian and (3=1 /UbT. The Lagrangian is a quadratic function 
of y and the coupling term between x and y is linear, so the functional integral over 
variable y can be performed explicitly by use of the Feynman- Vernon influence 
functional technique.^ Then the partition function is reduced to a single functional 
integral over x, i. e., 



Z = J T>x(t) exp 



S, 



2D r 



where the effective action is given by 

S%[x]=Sg[x]+Sg[x] 
r up 



5 eff L 



dr 



-Mx 2 
2 

r h/3 



Um(x) , 

Sg[x] = - jf dr J dr' [x(r) - x(r')f K(r r'). 



(11) 

(12) 
(13) 

(14) 



Thus the dynamics of the phase difference in a two-dimensional potential U{4>,4>s) 
can be mapped into one in an one-dimensional model. Note that due to the coupling 
between the junction and the LC circuit, the effective action S^g[x] contains a kind 
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of dissipation action <S£g[x] in a sense that the retarded (or nonlocal) effect exists. 
The nonlocal kernel K(t) in Eq. (fT4j) is defined by 

K(r) = ' -,\'„ - ^ - ( I > 



cosh 




\r\) 


sinh 


2 





In the next section, we will calculate the thermal and quantum escape rate from 
the effective action ■ 



4. Thermal Activation Process 

The TA escape rate well above the crossover temperature T co is given by ^ 



2D _ TT U l + + ^n7K) , , 

where is the quantum mechanical enhancement factor resulting form stable fluc- 
tuation modes, uj n is the Matsubara frequency, and 7(w n ) = (Cs/Cj)|a; ra |w£ C / 
w\ c ) is the Fourier transform of the memory-friction kernel. The potential barrier 
height is not changed even in the presence of the LC circuit, so V 2D = V^ 1D . Thus 
the coupling to the LC circuit modifies only the prefactor of Tta- 

For the non-adiabatic (u) p 3> ^lc) cases which can be applicable to actual 
experiment pl the quantum enhancement factor c 2 ^ (|17[) is approximately given by 



sinh 



C « — 7^ (18) 



In the case of Cs/Cj <C 1, the quantum enhancement factor coincides with 
the result without retardation effects, i.e., = sinh (h[3uj p /2) / sin (fr/3w p /2)PQ 
Therefore, the influence of the coupling to the LC circuit on the thermal activation 
process is quite weak, so the system behaves as one-dimensional systems well above 
the crossover temperature T co . In Sec. 6 we will numerically compare theoretical 
results with experimental data in the TA regime. 



5. Macroscopic Quantum Tunneling Process 

The MQT escape rate at zero temperature^ is given by 



'2 



r MQT = J im f - J ImlnZ (19) 
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By use of the bounce techniques, the MQT escape rate r^g T is perturbatively 
determined by 

2^ 



-B\y> — B v , 



(20) 



where £?id = 5*^ [xb]/7i = 36VQ LD /5fiw p and B rct = S^[xs]/fi are the bounce 
exponents, that are the value of the actions evaluated along the bounce trajectory 
xb{t). In the zero temperature and the non-adiabatic limit, the bounce action 
<S£ff [xb] is analytically given by 

2 



( u LC 



3 \ top 

Thus the total bounce exponent is given by 



LJpX l . 



157i 



(M + 8M)uj p x\, 



(21) 



(22) 



(23) 



where 

SM _ 5 Lj 
IF ~ 2 As' 

is the retardation correction to the mass M. Thus, due to the two-dimensional 
nature of the phase dynamics, the bounce exponent is increased with compared to 
one-dimensional cases. 

By substituting Eq. (|2"2"|) into Eq. (j2Tj)) . we finally get the zero-temperature MQT 
escape rate as 



t mqt( t = 0) = 



rkur> 



1 



6M\ 



36 Vr 



ID 



5 hui v 



1 



5M \ 

~m) 



(24) 



Therefore the coupling to the LC circuit effectively increases the barrier potential, 
i.e., Vq B — > V^ 1D (1 + SM/M). In contrast to the TA escape rate, the coupling to the 
LC circuit gives rise to reduce the MQT escape rate T^q T considerably. Therefore 
the two-dimensional nature has large influence on the MQT escape process at the 
low temperature regime. 



6. Comparison with Experiment 

In order to check the validity of the extended circuit 

mod ei™ 

for the TA and 

MQT escape process, we try to compare our result with the experimental data^ of 
the switching current distribution at the high-temperature TA and the low temper- 
ature MQT regimes. In these estimation, we will use Ic — 1.4 fiA, Ls = 1.7 nH and 
Cs = 1.6 pF which have b een directly determined from the TA and the microwave- 
assisted MQT experiments PE1 We have numerically calcul ated th e switching current 
distribution P(j) which is related to the escape rate T a tP^^ 



P( 7 ) = -r( 7 )exp 



(25) 
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Fig. 2. The temperature T dependence of the full width at half maximum cr of the switching 
current distribution P{"f). Both the one- (dotted black) and two-dimensional model (blue solid) 
give almost same result above the crossover temperature. The calculated a at T = K from the 
zero-temperature escape rate is shown by dashed-dotted (red) line. Experimental data of a (black 
circles) for a YBCO bi-epitaxial Josephson junction is also plotted. 



where v = \drj/dt\ is the sweep rate of the external bias current. In the ac- 
tual experiment,^ the temperature dependence of the full width at half maximum 
(HMFW) a of P{i) is measured. 

Firstly we will investigate the TA regime. In Fig. 2 we show the temperature de- 
pendence of cr in the TA escape regime (blue solid and black dotted lines). Both the 
one- and two-dimensional model give good agreement with the experimental data 
(black circles) well above the crossover temperature (T^ p ~ 0.04 K). Therefore, in 
the TA regime, the system can be served as an one-dimensional model without the 
LC circuit. 

In the MQT regime, the measured saturated- value of cr is found to be 11.9 nApl 
In order to numerically calculate cr, we need the information about Ic, Cj, Ls, and 
Cs- The values of Ic = 1.4 /iA, Ls = 1.7 nH and Cs — 1.6 pF have been directly 
determined from the experiments!^^ Therefore the only fitting parameter is Cj. 
From the numerical estimation of cr, we found that Cj = 0.22 pF gives good agree- 
ment with the experimental value of cr. The obtained value of Cj is consistent with 
the estimated value Cj « 0.16 pF based on the geometry of the junctionPSl There- 
fore, we can conclude that the extended circuit model can quantitatively explain 
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the MQT experimentP in the YBCO bi-epitaxial junction. 
7. Conclusions 

In the present work, the TA and the MQT escape process of the YBCO Josephson 
junction coupled to the LC circuit has been analyzed by taking into account the two- 
dimensional nature of the phase dynamics. Based on the Feynman- Vernon approach, 
the effective one-dimensional effective action is derived by integrating out the degree 
of freedom for the LC circuit. We found that the coupling to the LC circuit gives 
negligible reduction for TA escape rate. On the other hand, we also found that the 
MQT escape rate is considerably reduced due to the coupling between the junction 
and the LC circuit. These results are consistent with experimental result of a YBCO 
bi-epitaxial Josephson junction. 1121 In our model, we have assumed that the Josephson 
current-phase relation is given by Ij — Ic sin0. Therefore our theory can also be 
applicable to low-T c josephson junctions coupled to the artificial LC circuit. We 
expect that our prediction will be confirmed by use of not only high-T c but also the 
low-T c junctions experimentally. 
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